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Spectral and transport properties of the quasi two-dimensional adiabatic Su-Schrieffer-Heeger 
model are studied adjusting the parameters in order to model rubrene single-crystal field effect 
transistors with small but finite density of injected charge carriers. We show that, with increasing 
temperature T, the chemical potential moves into the tail of the density of states corresponding to 
localized states, but this is not enough to drive the system into an insulating state. The mobility 
along different crystallographic directions is calculated including vertex corrections which give rise to 
a transport lifetime one order of magnitude smaller than spectral lifetime of the states involved in the 
transport mechanism. With increasing temperature, the transport properties reach the Ioffe-Regel 
limit which is ascribed to less and less appreciable contribution of itinerant states to the conduction 
process. The model provides features of the mobility in close agreement with experiments: right 
order of magnitude, scaling as a power law T~ 7 , with 7 close or larger than two, and correct 
anisotropy ratio between different in-plane directions. Due to a realistic high dimensional model, 
the results are not biased by uncontrolled approximations. 

PACS numbers: 



I. INTRODUCTION 

In recent years, the interest in plastic electronics has 
grown considerably. In particular, organic field effect 
transistors (OFET) are the most common devices em- 
ployed to characterize the transport properties of the or- 
ganic semiconductors. Single- crystal FET constructed 
using ultrapure small molecule semiconductors have al- 
lowed to measure mobilities up to one order of magni- 
tude larger then those typical of thin film transistors. 1,2 
Among them, rubrene OFETs are much studied. 

Transport measurements from low to room tempera- 
ture in single crystal OFETs show a behaviour of the 
charge carrier mobility \i which can be defined "band-like" 
(/1 oc T~ 7 , with the exponent 7 close to two) similar to 
that observed in crystalline inorganic semiconductors^ 
The observation of the classical Hall effect is also ex- 
plained in terms of the presence of relatively free charge 
carriers, so that the charge conduction is incompatible 
with hopping transport between localized states/- 2 How- 
ever, the order of magnitude of mobility is much smaller 
than that of pure inorganic semiconductors, and the 
mean free path for the carriers has been theoretically 
estimated to be comparable with the molecular separa- 
tion with increasing temperature reaching the Ioffe-Regel 
limits 

Extended vs. localized features of charge carriers ap- 
pear also in spectroscopic observations. Angle resolved 
photoemission spectroscopy (ARPES) supports the ex- 
tended character of states^^ Indeed, photoemission ex- 
periments suggest that the quasi-particle energy disper- 
sion does exhibit a weak mass renormalization even if the 
width of the peaks increases significantly with tempera- 
ture. On the other hand, some spectroscopic probes, such 
as electron paramagnetic resonance (EPR), 6 ' 7 THz^ and 
modulated spectroscopy^ are in favor of states localized 



within few molecules. The theoretical explanation of the 
physical mechanism underlying the conduction proper- 
ties is still not fully understood. 10 

One of the main problems is to conciliate band-like 
with localized features of charge carriers^ In particu- 
lar, the nature of conduction at room temperature with 
mobilities close to the Ioffe-Regel limit remains contro- 
versial. Actually, in contrast with other polyacenes, in 
rubrene, to ascribe the presence of localized features to 
small polarons is not likely since the electron-phonon (el- 
ph) coupling is not large enough to justify the polaron 
formation. 10 

A model that is to some extent close to the 
Su-Schrieffer-Hceger (SSH) 12 one has been recently 
introduced. 11 It is a one-dimensional (ID) system where 
the effect of the electron-phonon coupling is reduced to 
a modulation of the transfer integral induced by low fre- 
quency intermolecular modes. This model applies to the 
most conductive crystal axis of high mobility systems, 
such as rubrene It has been studied by using a dy- 
namic approach where vibrational modes are treated as 
classical variables. This approach has been recently gen- 
eralized to two dimensions*^ Within this method, the 
computed mobility is in agreement with experimental re- 
sults. However, the role of dimensionality of the system is 
not clear: in fact, in the ID case, one has /i oc T~ 2 , while, 
in the 2D case, the decrease of the mobility with temper- 
ature is intermediate between \i oc T~ 2 and \jl oc T -1 . 
Moreover, this approach is not satisfying from different 
points of view. First, the dynamics of only one charge 
particle is studied neglecting completely the role of the 
chemical potential. Then, the dynamic disorder on the 
charge carrier due to the coupling with vibrational modes 
is included in a peculiar way: in fact, the electron is not 
coupled to classical gaussian fields with time dependent 
correlation functions which are determined by the con- 
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tact of the oscillators with an external bath^ In this 
case, the total system is composed only by one electron 
(or hole) and vibrational modes, therefore it is not clear 
if the corresponding coupled dynamics recovers the ther- 
mal equilibrium on long times and if the Einstein rela- 
tion can be properly used to determine the conductivity 
from the diffusivity. Finally, if the low frequency modes 
are not assumed much slower than the electron dynamics 
but occur in the same timescale, the quantum nature of 
phonons cannot be neglected. 

Recently, the transport properties of the ID SSH 
model have been analyzed within a different adiabatic 
approach^ Due to the geometry of OFET, in which the 
organic semiconductor is on the gate, one can assume 
that the oscillators are at the thermodynamic equilib- 
rium, therefore the problem is mapped onto that of a sin- 
gle quantum particle in a random potential (generalized 
Anderson problem). 17 However, in order to have finite 
results of conductivity, vertex corrections have been com- 
pletely neglected. Very recently, some of us have made a 
systematic study of this ID model including the vertex 
corrections?^ While finite frequency quantities are prop- 
erly calculated in this ID model, the inclusion of vertex 
corrections leads to a vanishing mobility unless an "ad- 
hoc" broadening of the energy eigenvalues is assumed. 

In this paper, we generalize a previous worbi^ in which 
some of usf have analyzed the properties of electrons cou- 
pled to local modes within the Holstein model by means 
of the adiabatic approximation^!^. We formulate an 
extension of the ID SSH model to the quasi 2D case 
since this is the relevant geometry for OFETs geometry. 
Moreover, the use of a high dimensional realistic model 
allows to include the physical anisotropy as an essential 
feature and to overcome the difficulties due to the lo- 
calization of all the states in reduced dimensions^ As 
discussed in Appendix, the proper scaling towards the 
thermodynamic limit is performed for the quasi 2D mod- 
els investigated in this paper. 

An important characteristic of our study is the small 
but finite carrier density with the introduction of the 
chemical potential. With increasing temperature, the 
chemical potential goes towards the tail of localized 
states. Therefore, all the experimentally quantities 
strongly dependent on the position of the chemical posi- 
tion will probe the features of localized states. The study 
of spectral properties clarify that the states that mainly 
contribute to the conduction process have low momen- 
tum and are not at the chemical potential. The mobility 
jj, is studied as a function of the el-ph coupling and parti- 
cle density. The behavior of mobility is in agreement with 
experiments. Not only the order of magnitude and the 
anisotropy ratio between different directions are right, 
but also the temperature dependence of fj, is correctly 
reproduced since it scales as a power law T~ 7 with 7 
close or larger than two. The inclusion of vertex correc- 
tions is relevant in the calculation, in particular to get 
a transport lifetime one order smaller than the spectral 
lifetime of the states involved in the transport mecha- 



nism. Therefore, high order corrections in two-particle 
correlation functions enhance the interaction effects on 
the properties of charge carriers. Actually, with increas- 
ing temperature, the Ioffe-Regel limit is reached since 
the contribution of itinerant states to the conduction be- 
comes less and less relevant. 

The paper is organized as follows. In section II, the 
model and approach are discussed in detail. The spectral 
and transport properties are discussed in section III and 
IV, respectively. Finally, section V contains conclusions 
and further discussions. 



II. MODEL AND CALCULATION METHOD 

The model studied in this paper include the anisotropy 
of small molecule organic semiconductors due to its high 
dimensionality. The following Hamiltonian defines the 
model: 

H = J2lku^+J2\Mu^+H eh (1) 
Ri Ri 

where mr ; is the classical oscillator displacement at the 
site in the position li^ the classical oscillator velocity 
at the site in the position Rj, M the oscillator mass, k 
the elastic constant. The whole lattice dynamic is due to 
a single effective phononic mode whose frequency ujq — 
yjKjM is assumed of the order of 5 — 6meVJ£ 

Since we are mostly interested in dc conductivity and 
low frequency spectral properties, we need to determine 
an effective Hamiltonian for the electron degrees of free- 
dom H e i in Eq. ([T] valid for low energy and particle den- 
sity. We start from the orthorhombic lattice of rubrene 
with two molecules per unit cellJ^ The dispersion law 
of the lowest among the two highest occupied molecular 
orbital (HOMO) derived bands is 

e (k) = — 2h a cos (k a a) — 2h b cos (k^b) — 2h c cos (k c c) 

fk a a + k b b\ (k a a-k b b\ 
—^hu+b cos I I - 2ha^b cos I i2) 

dropping the site energy diagonal terms. The quanti- 
ties a, 6, c represent the lattice parameter lengths along 
the three crystallographic vectors of the conventional cell. 
We assume a = 7.184, b = 14.443, c = 26.897^ We fol- 
low the experimental paper by Ding et al. in order to 
extract the transfer integrals. 4 The lowest HOMO de- 
rived band is measured through ARPES yielding the 
hopping parameters h a = 62.5meV, h b — l2.5meV, 

h a+b = ha-b = 57.5meV. The dispersion law at k c = 

22 

is drawn in figure (TT]). For small values of k a and k b , this 
graph strictly resembles a paraboloid. It is clear that an 
effective dispersion law 

e (k) = — 2t a cos (k a a) — 24 cos {k b b) 

is capable to fit well the dispersion law @ for small 
values of k a and k b and k c — 0. Fit estimates give 
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Figure 1: Dispersion law e (k) for rubrene at fc c = 0. 



t a = 118.6meV and t\, = 68.6meV. For t c there is no 
experimental measure but theoretical estimates seem to 
agree that it should be small compared with other direc- 
tions and it is likely owing to the large interplanar sep- 
aration of rubrene. In the following we assume t c much 
smaller than t a and As we will discuss later, vary- 
ing its value seems not too affect much the calculated 
mobilities. 

The electronic part H e i is then 



Ri,S 



a + -Or . 
R.+S 1 



(3) 



R, +8 



5 = a,b,c (4) 



In the last formula ap\ (aits) represents the creation (de- 
struction) of a charge carrier in the site R; . t ? is the bare 

effective transfer integral in the direction S, ag the el-ph 
parameter that controls the effect on the transfer inte- 
gral of the ion displacements in the direction 5. Once 
fixed a a , we impose ctb/a a = tb/t a and, in the same way, 

The dynamic contained in the term H e i is fully quan- 
tum, while we remind that the lattice obeys a classical 
dynamic. This assumption is justified by the adiabatic 
ratio huj/ta — 0.05. Therefore, the results discussed in 
this paper are valid starting from temperatures such that 
T > hu/kB, with ks Boltzmann constant. In the adia- 
batic limit the dimensionless quantity A 



(3) 



is the only relevant parameter to quantify the el-ph cou- 
pling strength. An ab-initio estimate 13 of Ao = 0.087 
has been given for rubrene and it has been substan- 
tially confirmed late n 22 ' 23 This represents our starting 
point, but a study of the mobility behaviour when A 
raises has been performed and is reported in the fol- 
lowing sections. In fact, if Ao is claimed to be a good 
value for the coupling in one dimension, keeping in 



mind the increase of the kinetic energy with dimension- 
ality, the equivalent A should be higher in three di- 
mensions. Simulations will focus on temperature sig- 
nificantly lower than t a /kB because this is the typical 
range investigated experimentally in OFETs. Moreover, 
in most OFETs the induced charge carrier concentration 
c = number of carriers / number of sites is rarely higher 
than 0.01 and this is the upper limit for the concentra- 
tions in our simulations. As shown in a previous pa- 
per dealing with a very similar but ID model ^ for this 
range of concentrations, the probability distribution of 
displacements is only slightly renormalized. Therefore, 
it is acceptable to refer to a free-oscillator (Gaussian) 
distribution for the displacements: 



P({un,}) 



_ ( 2w 



exp 



H 2 



<R, 



(6) 



where L = L a * Lb * L c is the finite size of the lattice 
and /3 = (&bT) _1 . We fix L c = 2 (i.e. we consider two 
crystalline layers) because in OFETs the effective channel 
of conductions cover very few planes^ Moreover, it is 
worth reminding that t c is very small. 

We use exact diagonalization to investigate both static 
and dynamical properties for the hamiltonian ([T]). At 
fixed configuration ({itR;}), one has to diagonalize the 
Hamiltonian yielding L eigenvalues E r , with r = 1, .., L. 
The eigenvector components Ui ir , with i = 1,..,L, are 
given through the unitary matrix U which diagonalizes 
the electronic problem. Each observable results from 
an average on lattice configurations. For example, the 
electron part of the partition function can be calcu- 
lated as (2 e i ({fR,})), where Z e i ({ur 4 }) is the quan- 
tum partition function of the electron subsystem at a 
given displacement configuration {hrJ. 25 Similarly, this 
would provide approximation-free evaluation of observ- 
ables such as mobility and spectral function in the semi- 
classical limits The real limitation is the accessible finite 
size of the lattice. The computing time for each diagonal- 
ization is of order L 3 and this constraints our analyses up 
to L a = Lb = 24. In order to reduce the finite size effect, 
we use periodic boundary conditions. We have checked 
that, for all the static and dynamic quantities studied 
in this paper, the thermodynamic limit is reached. The 
averages are performed by means of a Monte-Carlo pro- 
cedure. Actually, we generate a sequence of random num- 
bers distributed according to P ({ur ; }). For the systems 
investigated in this paper, to get a good accuracy even for 
dynamic quantities we perform the averages on a number 
of iterations up to ten thousands. 

Within the same framework, it is also possible to cal- 
culate the spectral functions, the density of states, and 
the conductivity. At fixed configuration, the density of 
states is D(w) ({ur ; }): 



(7) 



where Ar^r^w) ({urJ) represents the diagonal term of 
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the spectral function 

A Ri ,R» ({«rJ) = Krhrt (K - W) ■ (8) 



The spectral function in momentum representation 
Ak(uj) can be determined after performing the averages 
over the lattice configurations. In the case of conductiv- 
ity, for each configuration, we calculate the exact Kubo 
formula^ 



Re [<Tp, p (u) ({urJ)] 



7T (1 - 



■^p r (1 -p s 



r^s 



\(r\J p \s)\ 2 5{E s -E r + u) (9) 
where p — a,b 1 c, p r is the Fermi distribution 

1 



Pr 



1 + exp (j9(£7 r 



(10) 



corresponding to the exact eigenvalue E r at any chemical 
potential p p and (r |J p |s) is the matrix element of the 
current operator J p along the direction e p , defined as 



J p = ie 



E^( R *) 

Ri,S 



c n i c Ri+S' 



(11) 



with e electron charge. We notice that, in contrast with 
spectral properties, the temperature enters the calcula- 
tion not only through the displacement distribution, but 
also directly for each configuration through the Fermi dis- 
tributions p r . The numerical calculation of the conduc- 
tivity is able to include the vertex corrections discarded 
by previous approaches. 

Finally, the mobility is defined as 



lim 



Re [a p , p (u)} 



(12) 



For finite lattice sizes, the delta function appearing in 
spectral and transport properties has to be replaced with 
a Lorentzian, thus introducing a finite broadening rj 



6(E- 



1 



V 



u>) = lim 

r,^0+ 7T (E + UlY 



(13) 



Clearly, the limit procedure means that the value of rj 
should be kept as small as possible. This does not give 
any problem for quantities at finite frequency, as far as 
u) > rj. In the case of mobility, this procedure has to 
be correctly implemented since there is the limit u> — > 0. 
The correct expression is obtained performing the limits 
in the following order: u> t— > 0, rj and L i— > oo. These 
limits are not trivial, so that in the Appendix we have 
reported in detail the "calibration procedure" followed to 
fix the right rj for the calculation of the mobility. In 
contrast with the ID case, in our system the mobility 
dependence on the broadening rj exhibits a quasi-plateau 
with a maximum whose location gets closer to the zero 



as the size of the lattice increases. This makes possi- 
ble to set rj as the lowest energy scale of the system. 
In addition, it emerges quite clearly that, for the sizes 
considered in this work, the system is very close to the 
thermodynamic limit. This shows that the broadening 
•q is a parameter of the system which does not influence 
the mobility estimates. A detailed analysis has allowed 
us to obtain values for uj m i n and rj m i n that assure the de- 
termination of the correct value of the mobility for any 
temperature. For the mobilities presented in this paper, 
we used oj m i n = 10~ 3 < a , rj min of the order of 2 * 10~ 2 i Q . 

In the following, we will measure energies in units of t a . 
We will use units such that Boltzmann constant ks = 1 
and Planck constant h = 1 . We will consider the hopping 
parameters for rubrene derived in this section, therefore 
tb = 0.58. Along the c axis, we assume t c = 0.18. The 
small effects due to a change in the value of t c will be dis- 
cussed in the end of section about transport properties. 
Next spectral properties are analyzed. 



III. SPECTRAL PROPERTIES 

The study of spectral properties is important to indi- 
viduate the states that mainly contribute to the conduc- 
tion process. In Fig. [21 the Density Of States (DOS) 
is shown for different temperatures at A = 0.12 and 
c = 0.002. As shown in logarithmic scale, the DOS has 
a tail with a low energy exponential behavior. This re- 
gion corresponds to localized states^ We have checked 
analyzing the wave functions extracted from exact diago- 
nalizations that, actually, states with energies deep in the 
tail are strongly localized (one or two lattice parameters 
along the different directions as localization length) . On 
the other hand, beyond the shoulder (see Fig. [5] for Ti, 
T2 and T3), the itinerant nature of states is clearly ob- 
tained. This analysis gives rough indications for the mo- 
bility edge energy which can be located very close to the 
band edge for free electron E c (in our case, E c — —3.52). 

The states available in the tail increase with temper- 
ature. It is important to analyze the role played by the 
chemical potential ji p with varying the temperature. Ac- 
tually, ji p enters the energy tail and will penetrate into 
it with increasing temperature. At fixed particle density 
c = 0.002, for T = 0.12, one has ji p = -3.49, while, 
for T = 0.24, n p = -3.88 (see the inset of Fig. HJ}. 
One important point is that the quantity E c and the 
close mobility edge are always over ji p . Therefore, in 
the regime of low density relevant for OFET, the itiner- 
ant states are not at ji but at higher energies. We will 
show that these states are relevant for the conduction 
process. Therefore, the analysis of the properties of a 
high dimensional model points out that both localized 
and itinerant states are present in the system. This is a 
clear advantage of our work over previous studies in low 
dimensionalit y 11 ' 14 ^ 6 in which there are states that are 
more localized at very low energy and just less localized 
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Figure 2: The DOS as a function of the frequency (measured 
from the chemical potential /i p ) at A = 0.12 and c = 0.002 for 
different temperatures. E c is the band edge for free electrons. 
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Figure 3: The spectral function at k = as a function of the 
frequency at A = 0.12 and c = 0.002 for different tempera- 
tures. In the inset, the widths at half height of the spectral 
function is plotted as function of the temperature. 



close to the free electron edge. Apparently, in our sys- 
tem, with increasing temperature, more localized states 
become available close to and so the itinerant states 
become statistically less effective due to the behavior of 
the chemical potential. Eventually, the effect of pene- 
tration of /ip in the tail will overcome the enlargement 
of available itinerant states due to the Fermi statistics. 
We have checked that the trend of the chemical potential 
towards the energy region of the tail is enhanced with 
increasing the el-ph coupling A. 

The density of states can be calculated as the sum 
of the spectral functions over all the momenta k. 



We have checked that the spectral functions with low 
momentum are more peaked, while, with increasing k, 
they tend to broaden.— The tail in the DOS is due to a 
marked width of the high momentum spectral functions. 
In Fig. [3] we report the spectral function at k = for 
the same model parameters as the previous figure. The 
low damped states close to k = will keep their itiner- 
ant character, and they will be involved into the diffusive 
conduction process. Actually, we notice that the spectral 
function gives a negligible contribution to the weight of 
DOS at the chemical potential for all the temperatures. 
For example, at T — 0.20, the spectral weight is concen- 
trated in an energy region higher than that in which the 
chemical potential is located (fi p — —3.74 for c = 0.002). 
Actually, the spectral weight is in the region between 
fi p + 2T and fi p + 3T. Finally, we point out that, with 
increasing temperature, the peak position of the spec- 
tral function is only poorly renormalized in comparison 
with the bare one, in agreement with results of ID SSH 
models 

It is interesting to estimate the lifetime of the states 
at low momentum. The spectral functions at k = shift 
and broaden with increasing temperature due to the en- 
hanced role of the el-ph coupling. We have estimated the 
widths at half height in the inset of Fig. [3] At intermedi- 
ate temperatures, the width is of the order of 0.05t, that 
means of the order of 6 meV. Therefore, the lifetime of 
these states is of the order of 22 fs. In the next section, 
this lifetime will be compared with the transport lifetime 
derived from the mobility. 

The analysis of the spectral properties has clarified 
important features of the chemical potential and of the 
states participating at the transport process. Next, we 
investigate the "transition rules" between states involved 
into the conduction and analyze the mobility as a func- 
tion of the particle density and the el-ph coupling. 



IV. TRANSPORT PROPERTIES 

In this section, the focus will be on the Kubo formula 
for the conductivity given in Eq. (0). We will start from 
the analysis of the square modulus matrix elements of 
the current operator |(r | J a \ s)| 2 along the direction a be- 
tween exact eigenstates r and s. In Fig. [4j we report 
the contour plot of the current matrix elements in which 
the two axes are linked to the energies of the eigenstates. 
We consider two different temperatures T = 0.12 and 
T = 0.24. The dark color point towards high values of the 
matrix elements, the bright color to low ones. The dif- 
ference in the intensity between bright and dark regions 
correspond to square modulus matrix elements differing 
for about two orders of magnitude. The most important 
features emerging from the upper plot is that, in the en- 
ergy range important for low particle densities, there is 
a narrow region from about —3.4 to —3.3 in which the 
matrix elements give an appreciable contribution. More- 
over, states close in energy do not repel, but have sizable 
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Figure 4: Contour plot of square modulus matrix elements (in 
units of c2a2 /h 2 ) of the current operator along the direction a 
between exact eigenstates at A = 0.12 and T\ = 0.12 (upper 
panel), T2 = 0.24 (lower panel) . 



overlapping. These features are clearly ascribed to the 
role of itinerant states^ We notice that there are sev- 
eral gaps in the intensities between —3.3 and —2.6, then, 
very high matrix elements occur going toward the cen- 
ter of the electronic band. At higher temperature (Fig. 
31 lower panel), the contour map is broader, the gaps 
are almost completely removed and the absolute magni- 
tude of the matrix elements is lower. We point out the 
these "transition rules" are not strongly dependent on 
the particle density. We stress that, only for density c 
larger than 0.3 — 0.4, the matrix elements are sizeable, 
and the conduction states are far from the tail of the 
DOS, therefore one would get a "metallic state" typical 
of an inorganic compound. Indeed, both low density and 
el-ph interaction contribute to reduce the conducibilities 
of OFET in comparison with inorganic systems. 

The set of parameters used for the calculation of the 
mobility are discussed in Appendix A. Before calculating 
the mobility as a function of the temperature, we ana- 
lyze the mobility along the a direction as a function of 
the energy of the states which contribute to it. More pre- 
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Figure 5: The conducibility contribution along the a direction 
as function of the energies of the "ingoing states" at A = 0.12 
and c = 0.002 for different temperatures. The unit of the 
conducibility is fio where /Lto = ea 2 /h = 5.16cm 2 /(Vs). 



cisely, in Eq. we sum over the "outgoing states" s 
and analyze the results in terms of the "ingoing states" 
r. We again investigate the regime of low density rele- 
vant for OFET. In Fig. [5J we plot Ac/Aw as a function 
of w. At low temperature T = 0.12, the "initial" states 
relevant for the conduction are in the energy range where 
the current matrix elements are larger, i.e. from about 
u) = —3.4. Furthermore, they are in the range between 
l-ip+T and /i p +4T, so that they are itinerant. We, finally, 
note that the curve reproduces the same gaps as showed 
in Fig. [5] With increasing temperature, the itinerant 
states become more diffusive and distant from the chem- 
ical potential, and, within the energy range close to /i p , 
more and more localized states are present. This is the 
reason why the mobility gets lower at room temperature 
where it is nearly flat as function of the energy. 

In Fig. (J6j> , the mobility along the a direction as a 
function of the temperature is reported at fixed coupling 
A = 0.1 and for different concentrations c. The upper 
panel shows that the absolute magnitude of the mobil- 
ity substantially agrees with the experimental estimates 
being [i ~ 10cm 2 / (Vs) at room temperature, while, in 
the inset, where the scale is logarithmic, the mobility 
exhibits a quite linear dependence on the temperature 
that means a "band-like" behaviour with inverse power- 
law fi oc T~ 7 . The exponent 7 is evaluated as the slope 
of the straight lines drawn in the lower panel and the re- 
sults that we obtain from the fit are in the range 2 — 2.4, 
where the highest value is related to the lower concen- 
tration. This trend is in agreement with experimental 
measures that for rubrene establish 7 ~ 2 for tempera- 
tures T > 170 — 180K £ A feature of our model is that the 
mobility decreases for higher concentrations of carriers. 
This trend has been already found in the one-dimensional 
model, l = Actually, due to the "selection rules" discussed 
previously, increasing density can bring the system into 
an energy range where the current matrix elements get 
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Figure 6: Mobility along the a direction as a function of the 
temperature for different values of the concentration c. The 
coupling strength is fixed at A = 0.1. The inset of the lower 

panel reports the parameters of the fits via the function /i oc 
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Figure 7: Mobility along the a direction as a function of the 
temperature for different values of the coupling strength. The 
concentration c is fixed at c = 0.002. The inset reports the 
parameters of the fits via the function /i oc T 1 . 



lower. 

In Fig. ([7]) the mobility as a function of the temper- 
ature is reported at fixed concentration c and for differ- 
ent couplings starting from the value A = 0.087. We 
just notice two essential features emerging from Fig. ([7]). 
Quite obviously, the mobility decreases when the cou- 
pling strength is higher but, differently from the one- 
dimensional value which predicts the formation of bond 
polaron at A = 0.12^ in this case this does not occur 
even at higher values of A for which the "band-like" be- 
haviour still holds. To conclude this point, we notice that 
the effect of the coupling on 7 is less relevant then that 
of the particle density c. 

Starting from the mobility, we can determine the scat- 
tering time from the relation fi — eT tr /m. Since the 
mass is weakly renormalized from the el-ph interaction, 
one can assume m as the bare mass at k = 0. In the 
inset of Fig. [5J the scattering time as a function of the 
temperature is shown. We point out that is on the scale 
of the fs, so that it is one order of magnitude lower than 
the damping time of the states important for the spec- 
tral properties (on the scale of ten fs). Therefore, the 
transport processes amplify the effects of the el-ph inter- 
action and the vertex corrections introduced within our 
approach are fundamental to take into account this effect. 

From the scattering time, one can deduce the mean 
free path as l tr = v av T tr , where v av is an average velocity 
of the charge carriers. We estimate v av from the average 
kinetic energy as follows: 

Vav = J—(h)- (14) 

V m 

The operator (. . .) has to be intended both as a thermic 
average and as an average over the configurations {mrJ. 
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Figure 8: The mean free path (in units of the lattice param- 
eter a) as a function of the temperature at A = 0.12 and 
c = 0.002. The inset reports the scattering time as function 
of the temperature for the same model parameters of the main 
plot 



The quantity l tr is reported in Fig. [5] It is always on the 
scale of a few lattice parameters. The most important 
feature is its temperature behavior. As a consequence 
of the el-ph interaction effects, close to room tempera- 
ture, it becomes of the order of half lattice parameter a. 
This means that the Ioffe-Regel limit is reached^ The 
decrease of the mobility in the Ioffe-Regel limit is not due 
to a mass renormalization (dynamic and/or static) but 
is due to a reduction of the available itinerant states (the 
only ones able to transport current) with the tempera- 
ture. We remark that this result is due to the fundamen- 
tal role played by vertex corrections in the calculation of 
the mobility. 
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Another important quantity is the anisotropy of the 
transport properties along different in-plane directions. 
Up to now, we have discussed the mobility only along the 
a direction, since the corresponding quantity along the b 
axis can be roughly reproduced taking into account the 
anisotropy factors t /t a — 0.58. Actually, the mobility 
along b direction roughly scales as the square of the ratio 
tb/t a , therefore it is about 0.33 times the mobility along 
the a direction. This value is in good agreement with 
experimental results^ 

Finally, we briefly discuss the role of hopping parame- 
ter t c on the mobility features. Switching from t c = 0.18 
to t c = 0.10 produces a relative reduction of the mobility 
that lies between 10% and 15% depending on the con- 
centration of carriers c. Although the sensitivity of the 
mobility to t c is less relevant than other parameters, t c 
can't properly be treated as a "dead parameter". There 
may be two reasons for this effect: first reducing the in- 
terplanar coupling yields to a system more bidimensional 
and many theoretical predictions state that a purely bidi- 
mensional system should be an insulator^then, when t c 
becomes comparable with the broadening 77 (we keep 77 
in the range 0.1 — 0At c , see Appendix) the finiteness of 
the lattice appears partially invalidating the procedure. 
Probably, both effects give their contribution to the re- 
duction of the mobility. 



this result, it is important to include vertex corrections 
into the calculation of the mobility. The mobility as a 
function of temperature is in good agreement with ex- 
periments: indeed, the mobility has the right order of 
magnitude, it scales as a power law T -7 , with the 7 close 
or larger than two, and has the correct anisotropy ratio 
between different in-plane directions. The use of a real- 
istic quasi 2D model is fundamental for the analysis of 
this paper providing reliable results. 

The focus of this work has been on the effects of low 
frequency intermolecular modes on the transport proper- 
ties. While the mobility from 100K to 300A" is believed 
to depend almost entirely by this el-ph coupling, close 
to room temperature, the effects of other interactions, 
for example that due to local high frequency modes, can 
give a non negligible contribution.- - Actually, when the 
mean free path becomes of the order of the lattice param- 
eter as a result of the only intermolecular coupling, the 
local coupling can easily give rise to an hopping mecha- 
nism with an activated mobility. Rubrene is on the verge 
of this behavior, that can be accepted as a mechanism 
for other poliacenes. As a future work, it could be inter- 
esting to use the quasi 2D model discussed in this paper 
in combination with local el-ph coupling in order to bet- 
ter understand the conduction mechanism close to room 
temperature. Work in this direction is in progress. 



V. SUMMARY AND CONCLUSIONS 

Spectral and transport properties of the quasi two- 
dimensional adiabatic Su-Schrieffer-Heeger model have 
been discussed with reference to rubrene single-crystal 
field effect transistors. An important ingredient of our 
model is the small but finite carrier density, therefore the 
interesting behavior of the chemical potential as function 
of temperature has been investigated. Actually, with in- 
creasing temperature, the chemical potential always en- 
ters the tail of the density of states corresponding to lo- 
calized states. Therefore, all the experimentally quanti- 
ties strongly dependent on the position of the chemical 
position will probe the features of localized states. The 
combined study of spectral and transport properties is 
fundamental to shed light on the intricate conduction 
process of these materials characterized by small particle 
density and intermediate el-ph coupling. From this anal- 
ysis it emerges that the states that mainly contribute to 
the conduction process are the itinerant ones which are 
in the energy range between \i v +T and fi p + AT. With 
increasing temperature, these states are affected by en- 
hanced interaction effects and are more distant from fi p , 
and, at the same time, more and more localized states 
become available in the energy range close to the chemi- 
cal potential. Therefore, close to room temperature, the 
transport properties reach the Ioffe-Regel limit. More- 
over, the transport lifetime is almost one order of mag- 
nitude smaller than the spectral lifetime of the states 
involved in the transport mechanism. In order to get 



Appendix A: "Calibration procedure" 

In this Appendix, we report in detail the "calibration 
procedure" followed to fix the right value of the broad- 
ening 77 for the calculation of the mobility. According to 
Eqs. (|12p and (fT3|) . the mobility is obtained performing 
the two limits u> — > + and 77 — > + together with ther- 
modynamic limit L — > 00. It is clear that the limit on the 
size and on 77 has to be done at the same time because, 
actually they are similar. Moreover, it turns out that, for 
any fixed choice of the couple (77, L) and for any temper- 
ature T, the mobility depends slightly on the frequency 
to. There is a threshold under which the mobility is not 
sensitive to further decreasing of oj that means that we 
have reached the minimum frequency that our finite sys- 
tem can resolve. Practically, under the value ojq = 10~ 3 
the mobility shows a plateau. We have fixed this value 
for all the calculations of mobility. 

In Fig. ©, for three different values of temperature, 
we have plotted the mobility calculated via the Kubo 
formula ([9]) to analyze how it depends on 77 and the size 
L. The strength of coupling is fixed at A = 0.1 and the 
carrier concentration is c = 0.002. At any temperature, 
the series of data gets closer with increasing the size of 
the lattice. Eventually, for the two biggest sizes, they 
are not distinguishable for most of the range in which 77 
varies. This allows us to state that, for the maximum 
size that we can treat Lmax = L a * L * L c = 24 * 24 * 2, 
the system actually reaches the thermodynamic limit or 
at least gets very close to it. For example, in the worst 
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Figure 9: Dependence of the conducibility on r] for different 
sizes L of the lattice. Upper panel T = 0.12. Middle panel 
T = 0.18. Lower panel T = 0.24. The concentration of 
carriers is c = 0.002 for all the temperatures. The unit of the 
conducibility is fio where fio = ea 2 /h — 5.16cm 2 /(Vs). 



case of low temperature (T = 0.12 equivalent to 165K) 
the two upper curves split only for r) < 0.02t c and it is 
clear that the value tjmax for which the highest mobil- 
ity is reached gets closer to zero as the size rises. For 
all the regimes of temperature, the mobility exhibits a 
quasi-plateau (less definite at low temperatures) and the 
value of r/threshoid under which the mobility collapses is 
reduced by increasing the size. This trend can be easily 
explained if one notices that, when the size increases, the 
separation between two next eigenvalues reduces and the 
same happens to the width of the delta function required 
to couple a proper number of close eigenvalues. Another 
feature that is worth to be noticed is that the curves of 
different sizes merge together in better way at high tem- 
peratures then at low ones but, at the same time, for high 
temperatures T)max is higher. These arguments lead us 
to state that, for our model, a correct thermodynamic 
limit can be recovered in spite of the analogue ID model 
in which the mobility falls down under a value of 77 that 
in the limit of infinite size converges on a finite fjJ£. It 
seems that in ID the broadening 77 is not just a compu- 
tational need that virtually disappears in the infinite size 
limit but it's a real missing energy scale without which a 
finite mobility can't be obtained. Operatively, for each of 
the three temperature that appear in Fig. ©, we have 
determined the value rj ca i (T) for which the mobility has 
a maximum. Then, at any intermediate temperature T, 
the r] ca i (T) is obtained from a quadratic interpolation 
based on rf ca i (Tj). This procedure allows to calculate a 
mobility very close to that one reached in the thermody- 
namic limit. 
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